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Number  of  degrees  of  freedom. 

Unit  vector  parallel  to  L^. 

Unit  vectors  fixed  in  R.  (k=l,2,3) 

Number  of  cable  segments. 

A  typical  point  on  . 
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INTRODUCTION 


In  a  series  of  recent  papers  [1,2,3]*,  we  presented  a  method  for 
modelling  cable  dynamics.  The  method  is  based  upon  previously  developed 
general  procedures  for  finite  segment  modelling  of  multibody  systems 
[4,5].  In  this  report  we  present  results  obtained  by  applying  the  method 
in  a  series  of  submerged  cable  configurations.  The  results  are  com¬ 
pared  with  experimental  data  reported  by  Palo  at  the  Navy's  Civil 
Engineering  Laboratory  [6],  and  with  numerical  data  obtained  from 
SEADYN  [7] — a  finite  element  cable  computer  code. 

A  brief  review  of  the  method  itself  is  presented  in  the  first 
section  of  the  report. 


*Numbers  in  brackets  refer  to  References  at  the  end  of  the  report. 


I.  FINITE  SEGMENT  MODELLING 


We  model  a  cable  by  a  series  of  ball-and-socket  connected  rigid  links,  cr 
segments,  as  depicted  in  Figure  1.  The  dimensions  and  physical  para¬ 
meters  of  the  segments  are  arbitrary.  Also,  we  let  the  segments  be 
subjected  to  general  force  fields.  This  allows  us  to  simulate  an 
arbitrary  fluid  environment,  as  well  as  gravity  and  buoyancy  forces. 

Finally,  cable  stiffness  or  elasticity,  if  it  is  significant,  can  be  modelled 

by  springs  and  dampers  between  the  segments. 

.  «.  *  * 

We  describe  the  orientation  and  configuration  of  the  system  using 
the  relative  angles  between  the  segments.  We  assume  that  in  a  specific 
maneuver  that  one  end  of  the  cable  may  be  attached  to  a  towed  body  and 
that  the  motion  of  the  other  end  is  specified  (for  example,  it  may  be 
fixed).  We  assume  that  the  initial  configuration  of  the  system  is 
known.  The  objective  of  the  analysis  is  then  to  define  the  subsequent 
configuration  of  the  system. 


Figure  1.  A  Finite  Segment  Model  of  a  Cable. 
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Using  procedures  developed  for  finite  segment  modelling  of  multi¬ 
body  systems  {4,5],  we  can  obtain  explicit  governing  dynamical  equations 
of  motion  for  .the  cable  model.  What  is  more*  we  car.  obtain  these  equa¬ 
tions  in  a  form  suitable  for  conversion  into  algorithms  for  numerical 
integration.  Specifically,  the  equations  may  be  written  in  the  form: 


a.  X  «  f„ 

iP  P  l 


Ci.»P  "  !»•••*  n) 


where  the  Xp(p*l,...,n)  represent  the  relative  orientation  angles  of  the 
cable  segments,  where  n  is  the  number  of  degrees  of  freedom.  (Regarding 
notation,  there  is  a  s urn  for  repeated  indices  over  the  range  of  the  index.) 
a^  and  are  given  by  the  expressions: 


a,  «  m,v,  +  I.,  .m.  .0).., 

iP  i  jpk  I**  jkh.  jph  j£k 


f.  *  -  Kv.-.v,  ,  X  +  I £ 

i  i  j  jik  jqk  q  jkh  j£h  jqk  q 

+  e  ,  u i.  to,  u .  „  I .  k  i.  )  (3) 

mnk  jqn  jsr  j£k  jmr  q  s  K  J 

\  , 

where  is  the  mass  of  the  jth  cable  segment  ,  N 

is  the  n.trrber  of  cable  segments;  Ijjcjl(k,h“l,2,3)  are  the  components  of 
I j »  the  inertia  dyadic  of  with  respect  to  its  mass  center  Gj ,  referred 
to  unit  vectors  n^  fixed  in  an  inertial  reference  frame  R;  and  Vj 

Co=*l *  -  -  •  *N ;  p-1,...  ,n;  k*l,2,3)  are  the  n^  components  of  the  partial 
angular  velocities  and  partial  velocities  of  Lj  and  18];  F^(£=i, . . . >n) 

are  the  generalized  active  forces  [8]  developed  from  the  forces  applied 
to  Lj;  and  e^^  is  the  permutation  symbol  [9]. 


Equations  (1)  form  a  set  of  n  simultaneous  nonlinear  ordinary  dif¬ 
ferential  equations  determining  the  orientation  angles  Xp.  Since  the 
coefficients  a£p  and  f^  are  functions  of  the  arrays  vj£k  and  w  and 
their  derivatives,  the  system  of  governing  equations  can  automatically 
be  generated  once  the  and  uij ^  arrays  and  their  derivatives  are 

determined.  Simple  algorithms  for  computing  and  and  theiv 

derivatives  have  been  developed.  They  are  recorded  in  References 
[4]  and  (5] . 
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II.  FLUID  AND  GRAVITY  FORCE  MODELLING 


The  applied  forces  contributing  to  the  generalized  active  forces  F£ 
consist  of  fluid  and  gravity  forces.  The  fluid  forces  may  be  represented 
as:  1)  Normal  dr-  %  '"orcea:  2)  Tangential  drag  forces:  3)  Added  mass 
forces;  and  4)  Buoyancy  forces.  Reference  [2]  conrsins  a  detailed 
analysis  of  each  of  these  forces.  Specifically,  it  is  shown  in  [2]  that 
the  at  plied  force  at  ?  typical  point  P  (See  Figure  2.)  of  a  segment  L^ 
may  be  written  as: 


whv  -e  Vjj  xa  -He  v  naal  component  of  the  fluid  velocity  relative  to  the 
cable  segmer •  .  P,  is  the  normal  component  of  the  fluid  acceleration 
relative  to  the  cable  segment  at  P ,  is  the  tangential  component  of 
the  fl" ' i  velocity  relative  to  the  cable  segment  at  P,  w  is  the  segment 
weight  per  unit  length,  and  bjj  is  the  buoyancy  force  per  unit  length. 


Figure  2.  A  Typical  Cable  Segment  L^ . 


The  coefficients  A,  B,  and  C  are: 


A  *  C^p(r/4)d2 

(5) 

B  -  Cjjp(d/2) 

(6) 

C  -  CTp(d/2) 

(7) 

where  p  is  the  fluid  mass  density  and  d  is  the  diameter  of  L^.  C^, 

and  Qy  are  coefficients  dependent  upon  the  Reynolds  number  of  the  fluid 
flow  relative  to  the  cable  segment.  They  are  usually  determined  exper~ 
imentally  and  reported  results  may  vary  slightly.  Webster  [7]  records 
them  as 


-  1.0 


(8) 


and 


V 


0.0  for  ReN  <  0.1 

0.45  +  5.93/(Rejj)°*33  for  0.1  <  R£N  <  400.00 


.27  for  400  <  R^  <  105 

0.3  for  R  „  >  10s 
eN 


(9> 


1.88/(ReT) 


0.74 


for  0.1  <  Rey  £  100.55 


(10) 


0.062  for  ReT  >  100.55 


where  the  Reynolds  numbers  R  „  and  R  „  are  defined  as: 

en  ei. 


ReN  mpA\trJ/v 

(11) 

and 

ReT  “  pdlYTl/^ 

(12) 

where  y  is  the  viscoscity  of  the  fluid. 

The  weight  force  w  may  be  expressed  as: 

v  =  -ccgk 

(13) 
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where  pc  is  the  cable  ir,ass  per  unit  .length,  g  is  the  gravity  constant, 
and  k  is  a’  vertical  (up)  unit  vector. 

The  buoyancy  force  at  P  may  be  expressed  as  [2] 


b^j  *  p  ga  ,nx(k  x  n) 


(14) 


where  a  is  the  cross  section  area  of  and  p  is  a  unit  vector  parallel 
to  Lj.  (See  Figure  2.)  (Note  that  since  the  cable  segment  is  a  model 
of  a  portion  of  a  continuous  cable,  the  ends  of  the  segment  are  not 
exposed  to  the  fluid.  That  is,  the  segment  ends  are  "shielded"  from  the 
fluid  by  the' adjoining  cable  segments.  This  results  in  bN  being  normal 

If  the  set  of  applied  forces  at  all  points  along  are  replaced  by 
a  single  force  F\.  passing  through  Gj  together  with  a  couple  with  torque 
Tj,  then  F^  and  Tj  may  be  expressed  as: 


h  -  *M  +  +  ?T  +  S  +  h 

and 

h  m  h. +  h  . 


(15) 

(16) 


where  ?».  ?T*  W,  and  are  due  to  the  added  mass,  normal  drag, 
tangential  drag,  gravity  and  buoyancy  forces  respectively,  end  Tj^ 
and  Tjj  are  torques  due  to  the  added  mass  and  normal  drag  forces. 
These  forces  and  torques  may  be  expressed  as  [2]  : 


fM  =  AL?GN  (17) 

Fn  =  B(wx  n(l/3c)(X3/2  |_L/2  -  X3/2|l/2)  +  t(b/2c)«  x  n 

+  w[(Tfi>\/2  +  (£fef>HU/2i  ‘  1  <!8> 

+  Aac-lA  fiXn  + 

\  8c3/z  J  \A|_l/2-  (lA/2)  +  (b/2 
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and 


?T  -  “-I YgtIYct 

(19) 

W  = “vL  =  -p^gLk 

(20) 

-N  *  -NL  3  PgaL  5  x0s  x  S) 

(21) 

T„  *  -A(L3/12)  [a  -  (a  •  n)n  +  (m  •  n)n  x  u>] 

(22) 

IN  -  8S  *  Y^Ul/ScJX3'4  lL/2  -  (1/3c)X3/2  |.l/2  - 

b(cL  +  b)  vh\ 

8c  2  X  1 

=  Bn  x  yCTl(l/3c)X  -  lL/2  -  U/dcM  i_L/2  -  A  »L/2 

+  b(-cL+b)  ft  .  bte : -bi,loi  /*%2^W(b/2A  \ 

•  8c2  ~L/2  16cS"^  Vs  |  ,  ~  -(LcV)  +  <b/2c!5)/ 


_____  -L/2 

+  B[(w  •  e)n  -  u»K(l/4c)  [(L/2)  -  (5b/6c)  ]X3/2  { 

L/ 2 

+  (l/4c)[(L/2)  +  (5b/6c)]X3/^  [_L/2 
+  [5b2  -  4ac)  /  16c2  ]  *\/2  -  ^  *\l/2 

Vs I-  /9  +  (Lc’S/2)  +  (b/2cJs) 


(23) 


•  feir)  “-Sr 


L/2 

L/2 


-  (Lc^/2)  +  (b/2ci5)y 


where  L  is  the  length  of  L4  ,  <i>  and  a  are  the  angular  velocity  and  angu¬ 
lar  acceleration  of  L. ,  and  v_T  and  a_.  are  the 

j  -WJ  ~GN 

components  of  the  fluid  velocity  and  acceleration  normal  to  Lj  at  G,. 

Vgjj  and  a^  may  be  expressed  as : 

■  Vg  -  <Vc  •  B)5  (24) 

and 

~GN  =  ^W/G  “  (?W/G  #-S)S  (25^ 

where  Vg  and  §y/Q  ate  the  velocity  and  acceleration  of  the  fluid 
relative  to  G, .  Finally,  the  coefficients  a,  b,  and  c,  and  X  are: 

J  » 


R  -  V  •  V 

a  -GN  ~GN 

(26) 

b  =  ~2\T  *  u  x  n 

~GN 

(27) 
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and 


c  •  (y  x  r)  *  (y  x  n) 


(28) 
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III.  SUBMERGED  CABLE  DYNAMICS 


♦ 

Analytical  validation  of  the  finite  segment  model  (UCIN-CABLE) 
without  fluid  forces,  has  been  obtained  by  comparing  results  predicted 
by  the  model  with  results  obtained  using  other  methods.  ‘Reference  [3] 
records  details  of  this  validation.  To  obtain  additional  validation, 
and  to  develop  further  applications,  a  series  of  submerged  cable  man¬ 
euvers  were  simulated.  These  included  buoy  relaxations  and  achor  drops 
in  water  depicted  in  Figure  3. 

In  each  maneuver,  the  cable  length  was  12  in.  (1.83  m).  Silicon 
rubber  and  nylon  cables  were  simulated.  The  silicon  rubber  cable  had 
a  diameter  of  0.163  in.  (4.14  mm)  and  the  nylon  cable  had  a  diameter  of 
0.1  in.  (2.54  mm).  The  weight  densities  of  the  silicon  rubber  cables 
were  0.137  lb/in  (24.0  N/m)  and  0.108  Ib/in  (18.  92  N/m)  in  air.  (The 
heavier  cable  included  a  wire  conductor.)  The  weight  density  of  the 
nylon  cable  was  0.0448  lb/in  (7.85  N/a)  in  air.  The  buoys  and 
anchors  were  2  in.  (50.8  mm)  diameter  spheres  weighing  0.025  lb  (0.111  N) 
and  0.246  lb  (1.09  N)  in  air  respectively. 


Figure  3. 


Cable  Maneuvers. 


Six  tests  were  conducted.  The  first  four  were  buoy  relaxations  as 
in  Figure  3a.  The  last  two  were  anchor  drops  as  in  Figure  3b.  Twelve 
identical  cable-  segments  together  with  a  sphere  representing  the  buoy 
or  anchor,  were  used  in  the  model,  for  each  test.  The  fluid  mass  density 
and  viscoscity  were  given  the  values  1.9856  slug/ft3  (1024.16  kg/m3)  and 
3.516  x  10"5  slug/ft. sec.  (1.684  x  10“3  kg/m  sec)  to  simulate  seawater. 


The  tests  were  designed  to  simulate  experimental  tests  conducted  at 
the  Civil  Engineering  Laboratory  in  Port  Hueneme,  CA  as  recorded  by 
Palo  [6].  Table  I  provides  a  description  of  the  tests.  Figures  4.  to 
9.  show  the  cable  configuration  at  various  times,  together  with  compar¬ 
isons  with  experimental  results  and  finite  element  results  (SEADYN)  for 
the  respective  tests.  Figures  10,  to  15.  show  analogous  results  for 
the  fixed  end  tension.  Finally,  Figures  16.  to  21.  show  buoy  and  anchor  velocities. 


Test 

Number 

Type 

Initial  Position 
of  Buoy  or  Anchor 
(See  Figure  3.) 

Cable 

Material 

1 

Buoy  Relaxation 

x*51  in.,  z “-24  in. 

Silicon  Rubber 

2 

Buoy  Relaxation 

x«51  in.,  2 *-33  in. 

Silicon  Rubber 
with  Wire  Core 

3 

Buoy  Relaxation 

x*51  in.,  2 *-33  in. 

Nylon 

4 

Buoy  Relaxation 

x*51  in.,  z=16.6  in. 

Silicon  Rubber 
with  Wire  Core 

5 

Anchor  Drop 

x*  66  in. ,  2  “  0 

Silicon  Rubber 

6 

Anchor  Drop 

x*  54  in. ,  z  *  0 

Silicon  Rubber 

TEST 


Test  6:  Configurations  for  Anchor  Drop  for  Rubber  Cable. 
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TEST  1 


Figure  16.  Test  1:  Buoy  Velocity  for  Rubber  Cable. 


TOTAL  VELOCITY  OF  BUOY  (IN/SEC) 


TEST  4 


Figure  19.  Test  4: 


Buoy  Velocity  for  Rubber  Cable 
with  Wire  Core. 


IV.  DISCUSSION  AND  CONCLUSIONS 


In  the  current  set  of  tests  the  cable  was  relatively  light. 

Hence,  the  viscous  forces  were  relatively  large  compared  with  the 
gravitational  and  inertia  forces.  The  tests  thus  provide  a  validation 
of  the  fluid  force  modelling  of  the  UCIN-CABLE  code.  (As  noted  earlier, 
the  inertia  force  validation  is  reported  in  Reference  [3].) 

The  results  show  that  it  is  not  only  possible,  but  it  is  also 
practical  to  obtain  numerical  simulation  of  cable  dynamics  through 
finite  segment  modelling.  What  remains  is  a  validation  of  the 
modelling  for  three-dimensional  maneuvers. 
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